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The bispectrum of the cosmic microwave background (CMB) generated by a correlation between 
a time-dependent gravitational potential and the weak gravitational lensing effect provides a direct 
measurement of the influence of dark energy on CMB. This bispectrum is also known to yield the 
most important contamination of the so-called "local-form" primordial bispectrum, which can be 
used to rule out all single-field inflation models. In this paper, we reexamine the effect of non- 
linear matter clustering on this bispectrum. We compare three different approaches: the 3rd-order 
perturbation theory (3PT), and two empirical fitting formulae available in the literature, finding that 
detailed modeling of non-linearity appears to be not very important, as most of the signal-to- noise 
comes from the squeezed triangle, for which the correlation in the linear regime dominates. The 
expected signal-to-noise ratio for an experiment dominated by the cosmic variance up to i ma x = 1500 
is about 5, which is much smaller than the previous estimates including non-linearity, but agrees 
with the estimates based on the linear calculation. We find that the difference between the linear 
and non-linear predictions is undetectable, and does not alter the contamination of the local-form 
primordial non-Gaussianity. 



I. INTRODUCTION 

A time-dependent gravitational potential changes the 
temperature of the cosmic microwave background (CMB) 
as ST/T = 2 / dt (d^/dt) [1, 2], where * is a pertur- 
bation to the time-time component of the Friedmann- 
Lemaitre-Robertson- Walker metric. As ^ is constant 
during the matter-dominated era, a detection of this ef- 
fect directly shows that the universe is not completely 
matter-dominated, but has contributions from either spa- 
tial curvature or dark energy. Given the tight constraint 
on the spatial curvature we have from the current cos- 
mological data [3] , a detection of this effect is considered 
as the direct evidence for the effect of dark energy on the 
growth of structure [4] . 

The weak gravitational lensing effect caused by matter 
density fluctuations between us and the last scattering 
surface also changes the temperature of CMB, by shifting 
the observed directions of photons as T(n) — > T(n+d) = 
T(n) + d • VT(n) + . . . [5]. Here, the deflection angle d is 
given by d = 2 J dr r *~ r V i $>(r\r). As the same ^> enters 
in both effects, these two effects are correlated, yielding 
a non-zero 3-point correlation (bispectrum) in the CMB 
[6]. Therefore, this bispectrum can be used to probe the 
nature of dark energy [7]. 

The time-dependence of \& is caused by two effects: 
one is the linear growth, and the other is the non-linear 
growth. The former effect is caused by dark energy slow- 
ing down the growth of structure, leading to a decay of 
"J. The latter effect is caused by non-linear evolution of 



density fluctuations, leading to a growth of ty. Following 
the literature, we shall call the former the "integrated 
Sachs-Wolfe (ISW; [1]) effect," and the latter the "Rees- 
Sciama (RS; [2]) effect." In this paper, we shall calculate 
the bispectrum generated by the lensing-ISW correlation 
on large scales as well as the lensing-RS correlation on 
small scales. 

The main focus of this paper is the lensing-RS cor- 
relation. This correlation has been studied in the past 
[7-9] with different empirical methods for computing the 
non-linear matter power spectrum. In this paper, we sys- 
tematically compare two empirical methods used in the 
literature and the 3rd-order perturbation theory (3PT), 
which is applied to the lensing-RS bispectrum for the 
first time in this paper. We find that, while these dif- 
ferent methods yield somewhat different results for the 
lensing-RS bispectrum, the differences are too small to 
detect or affect our interpretation of the data. 

In Sec. II, we review the lensing-RS bispectrum. In 
Sec. Ill, we compare three methods for computing the 
non-linear matter power spectrum. In Sec. IV, we com- 
pare the lensing-RS cross-correlation power spectrum, 
Q(l), computed with different non-linear matter power 
spectra. In Sec. V, we calculate the expected signal-to- 
noise ratio of the lensing-RS bispectrum and the \ 2 dif- 
ferences between the linear model and various non-linear 
models. In Sec. VI, we study the effects of non-linearity 
on the contamination of the local-form primordial non- 
Gaussianity parameter, /nl- We conclude in Sec VII. 

Throughout this paper, we shall use the cosmologi- 
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cal parameters given by the WMAP 5-year best-fit pa- 
rameters (WMAP+BAO+P ML; [10]): Cl M = 0.277, 
n A = 0.723, h = 0.702, n s = 0.962, and a a = 0.817. 



II. LENSING-RS BISPECTRUM 

Let us use spherical harmonics to expand the observed 
temperature anisotropy, ST(h)/T = ^a; m Y/ m (n), as 
well as the "lensing potential," 6(n) = J2i m ©Zm^im(n), 
defined by d = VO. Then, the CMB bispectrum gener- 
ated by the lensing- RS correlation is given by [6, 11] 



hhl 3 



(fliimi Ct/ 2 m 2 dl 3 m 3 ) 



->mini2m3 



x C^(er 3m3 a^ 3 )+ Sperm. 



(1) 



Here, Cf is the primary CMB power spectrum without 
lensing, and ajf^ are the spherical harmonics coefficients 
of the ISW (or RS) effect. 

However, Eq. (1) is the leading-order contribution 
which is accurate only to ~ 10% level at I > 2000. 
Lewis, Challinor and Hanson [12] have shown that the 
sub-lcading-order correction can be incorporated by sim- 
ply replacing Cf above with the lensed CMB power spec- 
trum, Cf. 

h(h + 1) - h{h + 1) + h(h + 1) 
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We shall use this formula for computing the lensing-ISW 
(RS) bispectrum. Note that this prescription of replacing 
Cf with Ci is accurate only for the bispectrum in the 
squeezed configuration where one of the wavenumbers, 
say, l\, is much smaller than the other two, i.e., l\ -C 1% w 
1%. This is sufficient for our purpose because the signal- 
to-noise ratio of the lensing-ISW (or RS) bispectrum is 
dominated by the squeezed configuration. 
The Gaunt integral, G™f™* m3 , is defined as 
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Assuming statistical isotropy of the universe, rotational 
invariance implies that one can average over orientation 
of triangles (i.e., m's) to obtain the angle-averaged bis- 
pectrum [13]: 
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Here, the cross-power spectrum of the lensing potential 
and the ISW (or RS) effect, Q(l) = (@* 
by [6, 7] 

Q(i) = (e,* 



ISW 

3 m 3 -l 3 m 3 

z * r(z*)-r{z) dPy{k,z) 



2 Jo ^ r(z*)r{z) 3 



, (5) 



k=l/r{z) 



where z* = 1090 is the redshift of the last scattering sur- 
face, Pp the power spectrum of the Newtonian potential: 

P*(M)=nW) (j^j P(k,z)(i + zf, (6) 

and P(fc) the power spectrum of matter density 
fluctuations, 5m- This result follows from the 
Poisson equation (in natural units): k^f(k,z) = 
-4TrGp M {z)a 2 (z)5 M {k, z) = -A-KGp M o5 M {k, z)(l + z) = 
-lVL M Hl5 M {k,z)(l + z). 

On large scales where the scale-invariant spectrum, 
Pqi oc 1/fc 3 , is still preserved, the cross-power spectrum 
goes as Q(l) oc l/l 3 . On the other hand, the primary 
power spectrum goes as Cf oc l/l 2 . On smaller scales, 
Q(l) falls even faster than l/l 3 . This implies that the 
bispectrum peaks at the "squeezed triangle," for which 
one of I's is much smaller than the other two (e.g., 
I3 <C h ~ I2 if we order multipoles such that I3 < l 2 < h), 
and the smallest I corresponds to I of Q(l). This obser- 
vation suggests that the signal would be dominated by 
Q(l) in the small I for which matter fluctuations can still 
be treated as linear perturbations, and thus the detailed 
modeling of non-linear fluctuations may not be necessary. 
We will confirm this observation in this paper. 

The remaining task is to calculate P{k,z), including 
non-linear matter clustering. 



III. NONLINEAR MATTER POWER 
SPECTRUM 

A. 3rd-order Perturbation Theory (3PT) 

Higher-order perturbation theory is a promising ap- 
proach for computing non-linear evolution of matter den- 
sity fluctuations [14]. This is especially true at high red- 
shifts (z > 1), where non-linearity is not too strong [15]. 
The lensing-RS correlation has been studied using the 
3rd-order perturbation theory (3PT) by [16], who found 
a reasonable agreement between the 3PT prediction and 
the data obtained from the A-body simulation. 1 



1 Due to a page limitation, Ref. [16] did not report on the details 
of the 3PT results in the published version; however, the details 
are reported in arXiv:071 1.1696. 
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The matter power spectrum including the next-to- 
leading order non- linear correction is given by [17] 



P(k, 



[D(z)] 2 P 11 (k) + [D(z)} 4 [2P 13 (k) + P 22 (k)} 1 (7) 



where D(z) is a suitably normalized linear growth fac- 
tor (which is proportional to the scale factor during the 
matter era), Pn(fc) is the linear power spectrum at an 
arbitrary initial time, Zi, at which D{zi) is normalized to 
unity, and P 2 2{k) and Piz(k) are given by 



P 22 (fc) = 2 



where 



d 3 q 
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The 3PT is an attractive approach, as it provides the 
exact calculation in the quasi linear regime where the 
perturbative expansion is still valid. This should be con- 
trasted with the empirical approaches described below: 
they are calibrated using numerical simulations with a 
specific set of cosmological parameters, and thus cannot 
be easily extended to other cosmological models, such as 
dynamical dark energy models. 

A disadvantage of the 3PT is that its validity is lim- 
ited to the quasi linear regime, and thus the result on 
very small scales cannot be trusted. One can check the 
validity of the 3PT calculation by comparing it to a di- 
rect numerical simulation [15, 16]; or, one can compare 
an empirical formula calibrated to a specific cosmological 
model, to the 3PT calculation using the same cosmolog- 
ical model. We shall adopt the latter approach in this 
paper. 



B. Empirical Models 

Empirical approaches, which are calibrated using N- 
body simulations, have an advantage that they can, 
in principle, describe the matter power spectrum in a 
highly non-linear regime where the perturbative expan- 
sion breaks down. We shall use one of the popular meth- 
ods, called the "halo model," [18] for checking the va- 
lidity of the 3PT for computing the lensing-RS power 
spectrum. We then compare these results with another 
empirical model [19] used by most of the previous work 
on the lensing-RS bispcctrum. 



1. Halo Model (HALOFIT) 

In the halo model, the matter power spectrum is de- 
composed into two pieces: one that arises from two- 
point correlations between dark matter particles resid- 
ing in two different dark matter halos (2-halo term), and 
another that arises from two-point correlations between 
dark matter particles residing in a single dark matter halo 
(1-halo term). The former contribution is given approxi- 
mately by the linear matter power spectrum, whereas the 
latter contribution is given by the density profile of dark 
matter halos. This splitting between the 2- and 1-halo 
terms is somewhat artificial, and thus the halo model ap- 
proach would not provide an accurate description of the 
non-linear matter power spectrum, unless it is calibrated 
by numerical simulations. 

One popular calibrated formula is due to Smith et 
al. [20], which will be called "HALOFIT." They model 
the power spectrum as P(fc) = Pq(^) + Pff(fc), where 
Pq(k) is the quasi-linear 2-halo term: 



P Q (k)=P 11 (k) 



[l + fc 3 P 11 (fc)/(2^)] /3 " 



exp [-f(y)] 



l + a n k^P ll {k)/{2^) 

(11) 

where f(y) = y/4 + y 2 /8 with y = k/k a , and a n , /3 n , 
and k a are free parameters which need to be determined 
from simulations (see Appendix C of [20]). The second 
term, P^(fc), is the 1-halo term: 



Pff(fc) = 



a n y 



3/i 



1 + n n /y + v n /y 2 1 + b n yh + [c„/ 3 y]3-7. ' 

(12) 

where fi n , u n , a n , b n , c„, j n , fi, f 2 , and / 3 are free pa- 
rameters which need to be determined from simulations. 
As one of the cosmological models for which these func- 
tions are calibrated is a ACDM model with CIm = 0.3, 
CI A = 0.7, h = 0.7, and as = 0.9, which is close to the pa- 
rameters we adopt in this paper, this model can be used 
to check the validity of the 3PT results in the non-linear 
regime. However, we remind the readers that HALOFIT 
is not guaranteed to provide accurate results for the cos- 
mological models that are not explored in Ref. [20] , such 
as dynamical dark energy or massive neutrino models. 
Ref. [9] also used HALOFIT for computing the lensing- 
RS bispectruni. 



2. HKLM Scaling Model (MA 99) 

Another empirical formula is based on the idea origi- 
nally put forward by HKLM [21] for the real-space two- 
point correlation function. This idea has been applied to 
the power spectrum by [22, 23]. Then, Ma et al. [19] 
have extended the calibration to include models with dy- 
namical dark energy. We shall use the formula of Ma et 
al, and call it "MA99." 

HKLM postulates that the non-linear correlation func- 
tion is a universal function of the linear correlation 
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function, once the length scale (or the wavenumber) is 
rescaled by the mass conservation (i.e., transformation 
from Lagrangian to Eulerian coordinates). The form of 
this universal function needs to be found from numerical 
simulations. Ma et al. [191 find 



k 3 P(k,z) 
2^ 



G 



fc 3 P 11 (fc0^)/(27T 2 ) 



k$Pu(k ,z) 



2tt 2 



(13) 

where j3 = 0.83, <Ts(z) is related to the present-day as 
by crs(z) — [D(z)/D(0)](Xs, and go is defined by go = 



lw 



1.3ku|-0.76 



g(0). Here, g(z) oc (l+z)D(z), but it is nor- 
malized such that g(zi) = 1 during the matter-dominated 
era (e.g., Zi = 30). 2 The Lagrangian wavenumber, ko, is 
related to the Eulerian wavenumber, k, as 




0.001 0.010 0.100 1.000 

Wove number k [h Mpc -1 ] 



10.000 



ko = 



[l + k 3 P(k, Z )/(2TT 2 )} 1 / 3 ' 

The function, G(x), is given by 

1 + 0.02.T 4 + Cl x s /g(z) 3 



(14) 



G{x) = [1 + In (1 + 0.5x)] 



1 + C2X 



7.o 



(15) 



with ci = 1.08 x 10~ 4 and c 2 = 2.10 x 10~ 5 . 

How do we compute P(k, z) using this formula? 3 

1. Compute the linear power spectrum, Pn(fco, z), for 
a given Lagrangian wavenumber, fen. 

2. Compute G(x) in Eq. (13), and multiply it by 
k^P 11 (k ,z)/(2w 2 ) to obtain k 3 P(k, z)/(2tt 2 ). 

3. Compute the Eulerian wavenumber, k, using 
Eq. (14). 

4. Compute P(k,z) from k 3 P(k, z)/(2n 2 ) times 
(2^ 2 )/fc 3 . 



C. Comparing P(k, z) 

In Figure 1, we show non-linear power spectra at 
z = 0.1 and 1 computed from 3PT (solid line), MA99 
(dotted line), and HALOFIT (dashed line). We find 
that, somewhat surprisingly, HALOFIT and Ma et al. 
for z = 0.1 are fairly discrepant at k re 0.3 — 3 h Mpc -1 . 
(MA99 underestimates the power relative to HALOFIT.) 



2 For example, p(0) = 0.7646 for tt M = 0.277, fl A = 0.723 and 
w = — 1. 

3 While this formula has been used by most of the previous work 
on the lensing-RS bispectrum [7, 8], in all cases it has been im- 
plemented incorrectly. In the previous work, the authors used 
Eq. (13) with fcrj hi the argument of G(x) replaced by k. This 
is not the implementation proposed by the original paper [19], 
and goes against the original proposal made by HKLM. This 
observation gave an initial motivation for our work. 



FIG. 1. Comparison of non-linear power spectra computed 
from 3PT (solid line), MA99 (dotted line), and HALOFIT 
(dashed line). The linear power spectrum is shown by the 
dashed-dotted line. The upper and lower curves show P(k, z) 
for z = 0.1 and 1, respectively. 



In order to identify the origin of this discrepancy, we 
have also compared MA99 and HALOFIT with the for- 
mula by Peacock and Dodds (1996; PD96) [23] (not 
shown in Figure 1; in order to use their formula, it 
is necessary to use the smooth linear power spectrum 
without baryonic oscillations; thus, we used the smooth 
power spectrum given in [24]). We find that PD96 and 
HALOFIT agree well, which is consistent with the find- 
ing of [9]. However, MA99 and PD96, which are based on 
the same HKLM idea, differ significantly. This probably 
indicates that the difference already existed at the level 
of N-body simulations used by MA99 and PD96. Given 
that the latest HALOFIT formula has been shown to pro- 
vide excellent fits to a wide range of N-body simulations, 
we conclude that PD96 and HALOFIT are more accurate 
than MA99. This suggests that the previous work based 
on MA99 [7, 8] would require a reexamination. 



On the other hand, the 3PT results are close to 
HALOFIT, but lie slightly above it at k re 0.2 - 
1 h Mpc -1 . This is a known result: at a low redshift, 
the 3PT tends to overpredict the non-linear power spec- 
trum [15] . While this is an issue for the 1%-level precision 
cosmology using the galaxy power spectrum, the discrep- 
ancy at this level may not be so bad for the calculation 
of the lensing-RS correlation, as the statistical error on 
the expected total signal-to-noise of the measurement of 
the lensing-RS bispectrum is modest (S/N < 10). We 
thus take this as an encouraging sign and move on. 
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FIG. 2. Comparison of d\nP*(k, z)/dz for k = 0.1 ftMpc" 1 
as a function of z, computed from 3PT (solid line), MA99 
(dotted line), and HALOFIT (dashed line). 



IV. LENSING-RS CROSS-POWER SPECTRUM 

A. d\nP*{k,z)/dz 

The essential ingredient of the lensing-RS bispectrum 
is the lensing-RS cross-power spectrum, Q(l), defined by 
Eq. (5). In order to compute Q{1), we need derivatives of 
the potential power spectrum, Pj,, with respect to red- 
shifts, <9Pj(fc, z)/dz. This is related to derivatives of the 
density power spectrum, P(fc, z), as (see Eq. (6)) 

din Py(k,z) _ <91nP(fe,z) 2 

d~z ~ dz + (i + z y ( ' 

As P(fc, z) oc (1 + z)~ 2 for the linear matter power spec- 
trum during the matter-dominated era, dP^(k,z)/dz 
vanishes for this expected. When the universe is 

dominated by curvature or dark energy, the first term is 
still negative but becomes smaller than the second term, 
yielding dPy(k,z)/dz > 0. On the other hand, the 3PT 
result (Eq. (7)) shows that non-linear evolution gives a 
term in P(k,z) which goes as (1 + z)~ 4 , and thus one 
obtains non-zero <9Pj(fc, z)/dz even during the matter- 
dominated era. The sign is opposite: dP^(k,z)/dz < 0. 

We find that one needs to be quite careful about nu- 
merical accuracy when computing <91nPj(fc, z)/dz. A 
stable result can be obtained by the following method: 
compute lnPj,(fc,z) for various redshifts separated by 
5z = 10~ 2 , and then use a cubic spline interpolation to 
evaluate the derivative. For y(z) = lnPj(fc,z), 

t, , y(zhi) - y(zi ow ) 

y (z) = 

%hi ^low 

+ i[(3P 2 - l)y"{z hi ) - {3A 2 - l)y"(z low )](z hi - z low ), 

(17) 



FIG. 3. Same as Figure 2, but for k = 1 /iMpc -1 . 



where A = (z hi - z)/(z hi - z low ) and B = {z - z low ) / (z hi - 
zi ow ), and Zhi and z\ ovl denote the pre-computed values 
of redshifts that are closest to the chosen value of z. See 
Sec. 3.3 of [25]. This method gives a highly accurate 
9P*(fc, z)/dz compared to a simpler numerical differen- 
tiation such as y{z) = [y(z + Sz/2) — y(z — Sz/2)]/5z or 
y(z) = [y(z + Sz) — y(z)]/5z. We have verified this using 
the 3PT results: for 3PT, one can calculate the derivative 
exactly by differentiating Eq. (7) with respect to z: 

+2\D(z)] 2 \2P l3 (k) + P 22 (k)}} . (18) 

We find that the derivative from Eq. (18) and that from 
the cubic spline interpolation agree precisely. 4 

In Figure 2, we show 91nP*(fc, z)/dz for k = 
0.1 h Mpc -1 as a function of z, computed from 3PT (solid 
line), MA99 (dotted line), HALOFIT (dashed line), and 
the linear spectrum (dot-dashed line). At this wavenum- 
ber, they roughly agree with each other to within 5% 
at z > 0.01. Non-linear evolution of matter fluctua- 
tions makes d In Pj, (fc, z) /dz slightly smaller than the lin- 
ear prediction. At this wavenumber, the predicted non- 
linearity is the largest for 3PT, followed by HALOFIT 
and MA99, which is consistent with Figure 1. 

In Figure 3, we show 91nPj,(fc, z)/dz for k = 
1 h Mpc -1 . This regime is quite non-linear, and thus 
we see a clear change in the sign of d\nP<s,(k, z)jdz 
at a moderate redshift. (Recall that the linear evolu- 
tion due to dark energy gives a positive contribution to 



4 We suspect that any differences between our results presented 
in this paper and those presented in the literature [7-9] can 
be explained by either an incorrect implementation of MA99 
or an inaccurate computation of the derivative or both. A 
code for reproducing our results is available on http:/ /www.mpa- 
garching.mpg.de/~komatsu/CRL/ 
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FIG. 4. Absolute values of the lensing-RS cross-power spec- 
trum, |Q(i)|) as a function of multipoles, I, computed from 
3PT (solid line), MA99 (dotted line), and HALOFIT (dashed 
line). The linear power spectrum result, which does not show 
any change of the sign, is shown by the dashed-dotted line. 
Note that the sign of Q(l) is positive on large angular scales 
and negative on small angular scales. 



dlnP^ik, z)/dz, while the non-linear evolution gives a 
negative contribution to d\nP^(k, z)/dz.) The precise 
rcdshift at which the sign changes depends on models of 
non-linearity: it is z ~ 0.1 for 3PT while it is z ~ 0.3 for 
HALOFIT and MA99. 



B. Q(l) 

With dP^/dz computed, we now compute the lensing- 
RS power spectrum, Q(l), from Eq. (5). In Figure 4, we 
show \Q(l)\ computed from 3PT (solid line), MA99 (dot- 
ted line) and HALOFIT (dashed line). The sign change 
due to non-linearity is seen, and the multipole at which 
the sign changes depends on models of non-linearity. It 
is I ~ 640, 700, and 800 for 3PT, MA99, and HALOFIT, 
respectively. 



V. RESULTS 
A. Signal-to-noise ratio 

How well can we measure the lensing-RS bispectrum? 
The expected signal-to-noise ratio is given by [13] 



E 

2<l 1 l 2 l 3 <l a 



B 



C^C^Ci. 



E 

2<Zi<i 2 <i 3 <Z n 



B 



Aj 1 j 2 ! 3 Cz 1 Cz 2 Ci 3 



(19) 



FIG. 5. Expected signal-to-noise ratio of the lensing-RS bis- 
pectrum, S/N, as a function of the maximum multipole, Z max . 
All non-linear models as well as the linear model give similar 
results. 



where Aj^jj = 1 if all Vs arc different, A( 1 ; 2 ; 3 = 2 if two 
Vs are equal (isosceles configuration), and A; 1 ; 2 / 3 = 6 if 
all Vs are equal (equilateral configuration). 

This formula assumes that non-Gaussianity is weak, 
and the covariance matrix of the bispectrum can be ap- 
proximated by the Gaussian piece, C; 1 Cz 2 C; 3 . However, 
Lewis, Challinor and Hanson have shown that there is a 
non-negligible contribution from the no n- Gaussian signal 
generated by the lensing-ISW bispectrum to the covari- 
ance matrix [12]. We shall ignore this contribution for 
simplicity, as our primary goal here is to investigate how 
non-linear RS effect changes the signal-to-noise ratio rel- 
ative to the linear ISW effect. As a result, our signal- 
to-noise ratio is overestimated by 10% at l max = 1500 
and by > 40% at l max > 2000. As all the previous work 
except for [12] has also ignored this contribution to the 
covariance matrix, our results for the signal-to-noise ratio 
can be compared directly with those from the previous 
work. 

In Figure 5, we show S/N as a function of the maxi- 
mum multipole, / max . We find that all non-linear models 
as well as the linear model give similar results. This con- 
firms our earlier observation (see Sec. II) that the lensing- 
RS bispectrum peaks in the squeezed limit where the 
smallest multipole corresponds to I of Q(l), and thus most 
of the signal-to-noise comes from the region where Q(l) 
is well approximated by the linear lensing-ISW cross- 
correlation power spectrum. Our S/N estimate agrees 
well with that from the linear calculation of [26] . 

Note that our S/N for the linear model is about a 
factor of 2.7 smaller than that of [9]. 5 Their Eq. (15) 
suggests that they have not restricted the sum to li < 



5 For this comparison, we use Eq. (1) instead of Eq. (2) because 
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FIG. 6. The x 2 differences between the linear model and 
various non-linear models: 3PT (solid line), MA99 (dotted 
line), and HALOFIT (dashed line). 



FIG. 7. The contamination of /nl due to the lensing-ISW 
(linear effect only; solid line) and the lensing-RS (3PT; dashed 
line). Non-linearity does not affect the contamination of /nl- 
The solid and dashed lines are indistinguishable. 



h < h- which results in an overestimation of (S/N) 2 by 
a factor of 6, i.e., a factor 2.4 in S/N, which is enough 
to explain the difference. 

Our S/N for the non-linear model using MA99 is an 
order of magnitude smaller than that of [8] . 5 This is prob- 
ably due to a combination of their not restricting the 
sum to li < ?2 5~ ^3 and their overestimating the bis- 
pectrum with an incorrect implementation of MA99 (sec 
Sec. IIIB2). 

B. x 2 difference between linear and non-linear 
models 

Can we detect differences between the linear and non- 
linear models? In order to answer this question, we calcu- 
late the x 2 differences between the linear and non-linear 
models from [9] 

J2 ( B hhh ~ B hhh) / on \ 

where X and Y denote the names of models under con- 
sideration. For example, when we study the \ 2 differ- 
ence between 3PT and the linear model, X = 3PT and 
Y = LIN. 

In Figure 6, we show the x 2 differences between the 
linear model and various non-linear models: Xx-lin f° r 
X = 3PT, MA99, and HALOFIT. We find that, for all 
non-linear models, the x 2 differences are much smaller 
than unity, indicating that the differences are too small 



they have used the unlensed C; when calculating the lensing-RS 
bispectrum. 



to detect. We find similar values of x 2 differences among 
non-linear models. 

Our results do not agree with those of [9], who find 
y x-lin °f 01 'dcr unity for X = HALOFIT. We suspect 
that this is potentially due to (i) their not restricting the 
sum to l\ < I2 < ^3 (which would account for a factor of 
6), and (ii) a numerical accuracy of their evaluation of 
9P*(fc, z)/dz. As noted in Sec. IV A, a simple derivative 
such as dPy(k, z)/dz = [P*(fc, z + 5z) — Pm{k, z)]/5z can 
result in an inaccurate result, and a better method such 
as the cubic spline interpolation is needed for correctly 
calculating this derivative. We have confirmed this by us- 
ing the above simple derivative, finding that the results 
can vary significantly if such simpler numerical deriva- 
tives are used. 



VI. CONTAMINATION OF THE LOCAL-FORM 
PRIMORDIAL NON-GAUSSIANITY 

While the lensing-ISW and lensing-RS bispectra are 
useful for studying the nature of dark energy, they are 
also important because they yield the largest known con- 
tamination of the so-called "local form" primordial bis- 
pectrum [26-28]. The local- form bispectrum is particu- 
larly important, as a significant detection of the primor- 
dial bispectrum of this form would rule out all single-field 
inflation models regardless of the details of models [29] . 

The contamination of the local-form primordial bis- 
pectrum, parametrized by the /nl parameter, has been 
computed for the lensing-ISW bispectrum. How would 
non-linearity (lensing-RS) affect /nl? To answer this 
question, wc calculate the "bias in /nl," Le., a value of 
/nl which would be found if we fit the lensing-RS bispec- 
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trum to the local-form primordial bispcctrum template: 

^prim £jlens — RS 
1 ^ 2 ^ 3 Z 2 3 

<i<l a <l 3 A^C^Ct, 

o/nl = , m , (21) 

where Bf™™ is the local-form primordial bispectrum 
given in [11]. 

In Figure 7, we show the contamination of /nl for 
the lensing-ISW (solid) and the lensing-RS (computed 
with 3PT; dashed) bispectra. We find that they give 
similar results, and thus non-linearity does not affect the 
contamination of /nl- The values of S f^(l max ) that we 
find agree well with those from the linear calculation of 
[28]. 



VII. CONCLUSIONS 

The basic findings of this paper are simple: while 
non-linear matter clustering modifies the shape of the 
lensing-ISW bispectrum, differences between the linear 
prediction and non-linear predictions as well as differ- 
ences among non-linear predictions are too small to de- 
tect. Non-linearity does not affect the contamination of 
the local-form primordial bispectrum. This is because 
the lensing-ISW bispectrum peaks in the squeezed con- 



figuration in which the smallest multipole corresponds 
to the multipole of the lensing-ISW cross-correlation 
power spectrum, where the linear approximation is valid. 
Therefore, the linear calculation would be practically suf- 
ficient when interpreting the CMB data such as those 
from Planck. 

Nevertheless, if one wishes to improve upon the linear 
calculation, one should probably use the 3PT, as it offers 
a greater flexibility in terms of cosmological models for 
which the calculations are valid, as well as a straightfor- 
ward computation of dP^(k, z)/dz. 

Our results do not agree with the previous work study- 
ing the lensing-RS bispectrum [7-9] which found much 
greater effects of non-linear clustering on the lensing-RS 
bispectrum. We suspect that the discrepancy is due to a 
combination of an incorrect implementation of MA99, an 
inaccurate numerical evaluation of dP-^{k 1 z)/dz, and/or 
their not restricting the sum in (S/N) 2 to h < h < h- As 
a result, the expected signal-to-noise ratio of the lensing- 
RS bispectrum for a cosmic- variance-limited experiment 
is about 5 for l max = 1500, which is smaller than the pre- 
vious estimates [8, 9], but agrees well with the estimates 
based on the linear calculation [26-28] . 

We would like to thank J. Weller and D. Spcrgel for 
discussions. We would like to thank A. Lewis for pointing- 
out that we should use the lensed Ci in Eq. (2) as well 
as the importance of the non-Gaussian contribution to 
the covariance matrix given in Eq. (19). This work is 
supported in part by NSF grant PHY-0758153. 
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